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Abstract 

Within the framework of McLaughlin-Scott perturbation theory 
the equation-of-motion of the dislocation kink in the pinning potential 
is linearized, assuming the simultaneous influence of ac and dc forces. 
Based on the equations derived, the probability of kink depinning was 
studied. The dependencies of the depinning probability on dc and ac 
forces are analyzed in details. 

PACS: 61.72.Lk, 05.45.Yv, 05.40.-a 

1 Introduction 

In many elemental and compound semiconductors the motion of dislocations 
is controlled by Peierls mechanism, i.e. it occurs through formation and mi- 
gration of kink and anti-kink pairs. Formation of the dislocation kink implies 
that various dislocation segments are located in the neighboring valleys of 
the periodic relief (Peierls potential) of the lattice. Motion of a dislocation 
in Peierls relief is typically divided into three phases [I], [2]: formation of 
kink pairs, migration of the kink and anti-kink in opposite directions and 
annihilation of kinks with different signs moving towards each other. 

It is well known, that various point defects (impurity atoms, vacancies) 
may impede the motion of the kink along the dislocation, i.e. they can act 
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as pinning centers for the kink. In absence of lattice strain, the kink can be 
held in the potential of the pinning centre, while the influence of a stationary 
force can detach it with certain probability. Simultaneous effect of ac force 
can change (increase or decrease) the depinning probability. This effect was 
theoretically analyzed for the case, when a segment of the dislocation line 
overcomes the pinning centre [3J. 

The importance of the studied problem is mostly due to the issue of re- 
combination - enhanced motion of dislocations. In semiconductors the elec- 
tronic features of the dislocation kink can influence the dislocation dynamics, 
provided that the multiphonon capture of the charge carriers by the kink is 
the significant mechanism of their recombination. During the multiphonon 
capture the lattice fluctuations localized in proximity to the capture centre 
are stimulated. Subsequently, those fluctuations are dispersed through the 
entire crystal. Thus, if recombination of charge carriers involves the dislo- 
cation kinks and happens through the thermal capture mechanism [1], then 
the released energy excites additional fluctuations which may stimulate the 
motion of dislocation (similarly the photoplastic effect in semiconductors has 
been interpreted [S]). It is generally accepted that for dislocation motion by 
the Peierls mechanism, recombination basically contributes to the kink mi- 
gration along the dislocation [2J,[6]. Recombination - enhanced motion of 
dislocations leads to extension of stacking faults in junction regions of bipo- 
lar and optoelectronic semiconductor devices and eventually results in the 
device degradation. 

Thus, fluctuation-assisted depinning of the dislocation kink represents a 
single act that causes degradation of a semiconductor device. Therefore, the 
aim of our work is exploring stochastic dynamics of the pinned kink under 
the influence of stationary and alternating (harmonic) forces. During the 
degradation process the dislocation is influenced by constant force stemming 
from internal tensions of the lattice, and the lattice fluctuations resulting 
from thermal recombination of carriers act as the ac force. 

Kinks and anti-kinks on the dislocation are non-linear objects (described 
by the sin-Gordon (SG) soliton equation). Therefore, the study of fluctuation- 
assisted depinning of the kink is a complicated task. Simple approach to 
solving this problem is used in Ref. [7J, where dynamics of the kink is de- 
scribed by the equation of linear harmonic oscillator, which has the kink mass 
ascribed to it. This simplification allows to basically describe the process, 
however, for the purpose of further development of theory (to fully account 
for the dependence of detachment probability on the kink parameters), it is 
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necessary to correctly linearize the equation of motion. 

Linearization of the equation of motion of the kink is possible through the 
use of perturbation theory developed by McLaughlin and Scott [8|. Although 
this theory is based on some assumptions, it is actually the only reliable 
approach used in many studies for exploring the effect of external forces on 
the kink (see Ref. [H] and references therein). In this paper, we will use 
a modification of above-mentioned perturbation theory developed in ITU). 
When linearizing the equations of the kink motion in presence of etc force, 
this method enables to avoid complicated calculations used by the authors 



2 Kink dynamics under influence of ac force 

Let's consider a smooth kink on a dislocation (directed along x axis), i.e. 
the width w significantly exceeds the height a (Fig.l). The height a of the 
kink (equal to the inter-atomic distance) represents a period of Peierls relief. 
The kink is described by the transversal displacement (p of the dislocation, 
which is equal to zero when x —>■ — oo, and equal to a when x —>■ oo (for 
anti-kink the reverse is true, i.e. ip decreases from a to zero). The kink width 
corresponds to the range where < tp < a. The smallness of the parameter 
Tp/G is the criterion determining the smoothness of the kink, where Tp is 
the Peierls tension and G is the shift module. The kink width relates to its 
height as per following equation rjQ: 



Therefore, the inequality Tp/G << 1 denotes that the smooth kinks are 
being formed, as from ([T]) we get the following: w » a. In this case the 
string model of the dislocation is applicable for description of the kink. 

If the kink is located close to the pinning point, we can describe the 
interaction of the kink with the pinning centre (located at the point x = 
on the dislocation) by the following potential [8j,[10j: 



where parameter /z characterizes the interaction force. The multiplier 1 — 
cos(27rip/a) in equation (T2J) ensures that the pinning potential U becomes 



[8j. 




(1) 
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equal to zero outside of the kink region. 

As indicated in the Introduction, we are interested in the kink dynamics 
under influence of dc and ac forces. The constant force T stemming either 
from the internal pressure or the lattice external strain tends to move the 
kink away from the pinning centre. We will consider the ac force 7(i), which 
emulates the lattice fluctuations originating due to the recombination pro- 
cess, as harmonic process. Assuming that the energy of these forces, as well 
as dissipation due to friction are small compared to Peierls potential, we 
can consider the influence of external forces as small perturbation relative 
to the undisturbed system. If the periodic potential has the following form 
(commonly accepted for Peierls relief [Tj.|llj): 

Vp = *£*?(m, ( 3 ) 



27r V a J 

then the dynamics of a string in this potential in case of existence of dissi- 
pation can be described by disturbed SG equation [8]: 

na X . 2ir<p 

nuptt ~ XVxx = —— o sm h /, (4) 

2w z a 

where m and x are the mass and the linear tension of the dislocation (per 
unit length) respectively. The force / is the sum of all acting forces: 

/ = -r}ip t -r - "f(t) + fj,5(x)sm— -, (5) 

a 

where rj characterizes dislocation damping per unit length (rj > 0). For 
derivation of equation (TJJ it was taken into account that \ — Gb 2 /2 [TT], 
where b is the module of the Burgers vector (which has the order of a), s is 
the sound velocity in solids, and the Tp/G ratio is determined by equation 

©• 

If the system is not perturbed, i.e. / = 0, then the solution of equation 
(HD is the following 0,0,110]: 



2a 

<p = — arctan 

7T 



exp (±—g (v) (x — xq — vt)^j 

where Xq is the initial coordinate of the kink, v is the velocity of motion across 

the dislocation (v < s), and g(v) = 1/yl — {v/ s) 2 . Here, the positive sign 
corresponds to the kink and the negative sign - to the anti-kink. This solution 
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corresponds to translation invariance of the kink along the dislocation. The 
energy of SG system is presented by the following Hamiltonian [5].[TU]: 



V 2s 



w 



By substituting the solution of unperturbed SG equation into © and inte- 
grating over x, we will arrive at the Hamiltonian of the unperturbed system: 

» sc =^(i-w*)y 1/2 . m 

If / << arp/2, i.e. all forces are small compared to the amplitude of the 
Peierls potential, then in order to solve equation (j3J) we can use the pertur- 
bation technique developed in publications [5], [ID]. We assume that solution 
of equation (j3J) has the form of unperturbed SG equation, however the per- 
turbations lead to modulation of velocity, i.e. in this case v is a function of 
time t: 



2a 

ip = — arctan 

7T 



exp (±^(g(v)x-X(t)f) 



(8) 



where X (t) — J g (v (t')) v (f) dt' indicates the location of the kink centre 

o 

(assuming that x$ = 0). In absence of terms describing dissipation and ac 
force in equation (jSJ), we can write the Hamiltonian of the perturbed system 
the following way: 

H((p) = H SG + iP, (9) 



where 



oo 



H* = - [ (Tip-^-fiS(x)(l-co8^))dx (10) 



w J V 2tt 



-oo 



describes the energy of stationary perturbation. In this case H(<p) is station- 
ary in time: dH/dt = 0. When non-stationary forces are present (the friction 
force and ac external force), the equation of motion must be determined by 
the following equation 



/ (w 2 + j{t)<p)dx. 11 

dt W J-oo v ' 
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By inserting equation into (11 II) . we arrive at the equation of the momen- 
tum of SG system: 



dH SG 7T 7 

——i, = - / f(<p)<pdx, (12) 



-oo 



where cp = <p (g (v) x, X) is the solution of the perturbed SG equation (jlj) 
that have the time derivative (p ~ ipxg{v)v [10J. The equation ffl2l . along 
with the equation of time derivate of X 

X = g(v)v, (13) 

describes the motion of the kink. The Hamiltonian H SG in equation ffl2|) has 
the same form as the Hamiltonian of the unperturbed SG system ([6]). The 
string displacement <p is determined by (jSJ). By inserting equation (jSJ) into 
(IHD and integrating, we arrive at the expression for the soliton energy (which 
substitutes the stationary Hamiltonian (j7j)): 

H SG [p (g (v (t))x, X(t),v (t))) *> (i_ {v {t ) ft 2) - 1/2 . ( M ) 

Now, by inserting the expressions (TH1) and (jSJ) into ffl2|) and integrating, we 
get set of differential equations for variables v (t) and X (t) 

Mf t = ± [r + 7 (*)] (^))~ 3 - ^ (5 (^))" 2 - 

%(g(v)y 2 sech>(?£)unh(?£) , (15) 

where we accounted for x — ms 2 ; M = 2am/ 'nw is the effective mass of the 
kink. The set of equations (fl5j) determines the dynamics of the isolated kink 
(i.e. the one which does not interact with other kinks) conditioned by the 
general force /. 

Let's explore the motion of the kink in the vicinity of the pinning point. 
Equations (Jl5j) are typically analyzed in the phase plane , X). The singular 
point, where v — and X = 0, in presence of dc force corresponds to the kink 
equilibrium position [8]. When ac force is also in effect, we have to deter- 
mine the equilibrium position considering only stationary force, meanwhile 
the ac force causes additional shift from the equilibrium. In other words, the 
equilibrium position is determined as the kink coordinate averaged on the 
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time intervals larger than the period of the ac force. As it follows from equa- 
tion (1151) . the following conditions correspond to the position of equilibrium: 
v = and X = Xq, where Xq is the root of the following equation: 

f-*- 2 (^H(*H («•) 

For the purpose of analytical examination of the kink oscillations around the 
point Xq, we need to linearize the equation ffl5l) with respect to variable v, 
assuming that the motion of the kink is non-relativistic (v « s). In absence 
of ac force, the kink makes damping oscillations around the pinning point, 
and eventually reaches the rest [Hj. If the kink is also subject to the influence 
of periodic force 7 (t) = fo sin fit, where fo is the amplitude of periodic force 
and Q is its frequency, then we need to analyze the following equations: 



Mf t = ± [r + /o sin fit] - %v - % seek* (g) tanh (^) 

dX _ . ' \ / 



dX 
dt 



Further linearization of the equation of motion of the kink is possible in 
case of small oscillations, assuming that the displacement is small compared 
to the kink width (X — Xq << w). Then we can expand the last term 
in the right-hand side of the first equation (117p (describing the interaction 
of the kink centre with the pinning point) into Taylor's series around the 
equilibrium point Xq. Keeping only linear in X — Xq term and using the 
second expression ffTTj) . we obtain the equation of linear oscillations: 



X + 2XX + LO^(X-X ) = ±jtcosm, (18) 
where the parameter A = arj/7rwM is the effective damping coefficient, and 
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w 2 M 



4 {nX \ n _ L2 f7rXQ\^ ^ 2 f7rX 



sech —2 sech tanh 

V w ) V w J \ w 



(19) 



assumes the role of eigenfrequency of the system. 

As it follows from equation when T = the equilibrium position 

is Xq = 0, and from equation (TT9]) we derive Uq = Jfp. With increasing 
constant force, c^o reduces; the critical value of force r c = 4/x/ (3y/3u?J corre- 
sponds to the condition ujq = 0, which is determined from equations ( Tl6|) and 
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(TT^|) . In other words, when T > T c (then equation (1161) haven't roots), the 
kink overcomes the pinning potential and can be considered as a free kink. 

Note that the constant force does not enter explicitly into equation fjlSp . 
This force determines the equilibrium position Xq. Thus, equation ffl8|) de- 
scribes oscillations of the shifted coordinate of the kink centre under the 
influence of the harmonic force. In case of small viscosity ujq > A, we imme- 
diately arrive at the kink displacement 

X = X + s(t), (20) 

where the function s (t) = B cos (Qt + 5) describes harmonic displacements; 
B is the amplitude of the kink oscillations and 5 is the initial phase: 

B = = .tano 



M v /(^-n 2 ) 2 + 4A 2 fi 2 ' W-ul 



3 Fluctuation-assisted depinning of the kink 

We will use Langevin approach to describe fluctuations of the kink in the po- 
tential of the pinning centre, i.e. fluctuation force f(t) needs to be introduced 
into the right hand side of the equation [T3] : 

X + 2XX + u 2 (X - X ) = ± A cos Qt + (21) 

We consider the fluctuations as being white noise, i.e. the average value of 
the fluctuation force equals to zero, and the time correlation is determined 
by the delta-function: 

(/») = o, (22] 

(f(t)f(t>)) = 2D8(t-t>), { ] 

where D is the intensity of the Langevin source (it also represents the coeffi- 
cient of diffusion in velocity space). From now on, we will not consider quan- 
tum fluctuations (i.e. we assume that the process temperature is higher than 
the Debye temperature of the crystal T >> 9 D ); in this case D = Xk B T/M, 
where ks is the Boltzmann constant. 

Before proceeding to solution of the Langevin equation, let's determine 
the value of the critical shift of the kink centre, which corresponds to the 
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kink depinning. Obviously, a parabolic potential corresponds to the lin- 
earized equation (TTBT) . However, we have not yet determined the radius of 
that potential corresponding to the critical shift. In order to determine that 
radius, let's refer to first equation ffTTj) . We can determine the local potential 
of the interaction of kink and pinning centre U(X) by integrating over X the 
last term of this equation (characterizing the strength of interaction). Taking 
into account that this potential should become zero when X — > oo, we get 

[/(X) = -t sech 2 ( —) . (23) 



7T V W 

Unfortunately, the potential ( 1231) does not have the natural cutoff radius. 
Therefore, let's assume that some small value of the interaction potential U = 
—Uq corresponds to the critical shift value X cr . The kink can be considered 
as detached from the pinning centre if its energy state is separated from 
the zero level (corresponding to the free kink) by the energy Uq ~ ksT: in 
this case any single- quantum thermal fluctuation will lead to final exit from 
the potential well. Thus, critical shift can be determined by the following 
expression: 

Now we can proceed to calculation of the probability of the fluctuation- 
assisted shift of the kink beyond X cr value ( l24l) . The solution of equation 
( I2TT) is the sum of the deterministic shift (118]) and shift £ (t) determined by 
fluctuations: 

X = X + s{t)+£{t). (25) 



v w h 
X cr = — arcosh 



(24) 



By using the theory of random process overshoots [14J, we can determine 
the depinning probability of the kink affected by the periodic force as equal 
to the average (during oscillation period) speed of positive overshoots of the 
sum of stochastic and harmonic shifts beyond the value C = X cr — X : 

W = — J dt JdvvW{C-s(t);v-s (t)) (26) 



where W is the joint probability density for the coordinate and velocity of 
the random parameter. Since we analyze influence of the white noise f(t) on 
the system described by linear equation (1211) . then the response f )25|) should 
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be a normal process described by Gauss function: 



W(X,v) = (27ro-o-i) -1 exp 



(*-«( *)) 

2(X 2 
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exp 



2a\ 



21 



(27) 



where a and 01 are the dispersions of the coordinate and the velocity of 
the fluctuating parameter respectively, which are determined through the 
correlation function of the process. 

The spectral density of the process S (a?) = (X 2 (a?)) is determined from 
equations (12T1) and (1221: 



5 M 



D 



(cu 2 - c^ 2 ) 2 + 4A 2 c^ 



and for the correlation function K (t) = j- J S (u) e luJT du we find: 



K ^ = 2~M7 exp ' T '^ ' cos ^ ~ A ~ r 



A 



ul - A 2 



sin \/a;o — A 2 |t| 



(2f 



Dispersions of the coordinate and the velocity are determined by a = K (r = 0), 
a x = K"(t = 0): 

2 irksT 2 nksT 



By using expressions (l27j) and (|29|) . we can now perform integration over v 
in (I26I) . By denoting by \l/ the phase of deterministic shift = Qt + 5 and 
considering that integrand is periodic as function of ^ with the period 2tt, 
we will get 



fa 



-J^J $'(c-pcos^) X F X (-^;-;--(pgsin#) 2 )d*, 



(30) 



where $' (2) = ;^=exp (— 2 2 /2) is the derivative of the probability integral, 
1-F1 (a; b; z) is the degenerate hypergeometric function; p = B/a, q = £1/uJq, 
and c = C/a is the relative level. 

In order to explore the dependence of the kink depinning on the stationary 
force, we discuss first the case when ac force is absent. In the latter case, for 
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the probability of the kink depinning the following expression can be easily 
obtained: 

" o = S exp (-y)- (31) 

Dependence of v$ on constant force V is determined through dependence of 
the frequency u a and relative level c = \J ™^t (X cr — X ) on constant shift 

Xo (while the dependence of Xo on V is determined by the transcendental 
equation (fl6l) ). The analytical function vq (T) can be obtained in case when 
constant force T is significantly smaller than the critical value of T c (r << 
r c ). By denoting y = tanh (irX /w), from equation ( fl6l) we can get the 
cubical equation: 

or 

:i ' :l " (32) 

In the case T << T c it has the following solution 

2 / r 



V = 1= sin . 

By limiting to the first order of value T/T c , we will get the following depen- 
dence of Xo on the force: 

X = --artanh ^-j. 

By inserting this expression into ffl9|) . we will find the relationship between 
the eigenfrequency of the system and dc force: 

*-3£H (£)")■ 

Finally, using equation (133!) the depinning probability is obtained: 



V27ru; 2 M/ I 3 V 3f r 



exp 



2u / i /4r \ 2 \ / , / /i \V2 2 r \ 2 

arcosn — — — + 



7i 2 k B T\ 3 \3r c y J \ \7Tk B Tj 3V3T 



(34) 
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The dependence of ( 1511) is presented on Fig.2. Note that for a value of dc 
force T (r < T c ) which is determined by the equation 



H£ - sech* (^) tanh (?*~) = 0, 
2u \ w J \ w J 



the level C = X cr — X becomes zero. Then the depinning probability during 
the period 2ti/ujq approaches unity as it is shown by dashed curve on the 
graph. 

The dependence of probability fl30|) on the amplitude of ac force B = pa 
can also be obtained more explicitly in several limiting cases. Note that the 
condition of incomplete equilibrium q = Q/ujq << 1 should be fulfilled for 
the quasi-stationary processes under discussion [T5]. Therefore, depending 
on the force amplitude, two limiting cases can be realized: pq << 1 and 
pq » 1. 

a) The small amplitude of the ac force, pq « 1. In this case we can 
derive from equation (1301) an expression linking uq with the probability of 
detachment in absence of ac force u : 



The dependence of vq (p) for various values of the relative level c is depicted 
on Fig.3. In the case c < 1 , Uq is a monotonous decreasing function relative 
to signal p, i.e. if the level of C is less than the amplitude of noise a, then 
presence of external periodic force reduces the depinning probability. 

If c > 1, vn{p) increases reaching maximum when p = p m (c): p m = 
2y/l ~ c" 2 for pc < 1, and p m w c for pc > 1. Thus, when interaction 
strength \x ensures high values of critical shift X cr (and respectively, the level 
C exceeds the noise amplitude) and p ~ p m , the probability of detachment 
can significantly increase. When signal continues to grow (p > p m ), the 
probability reduces as uq ~ p^ 1 . 

b) The limit of large amplitudes, pq » 1. In this case the probable case 
is p > c. As it follows from ( |26l) and (1271) . uq — > Q/2tt, i.e. the number of 
overshoots during the period approaches unity. 

Note that the dependence of kink depinning probability on amplitude 
of ac force is analogous to the case of dislocation line (string) detachment 
from the oscillating stopper |3j. The latter is true for dislocation motions in 
crystals with low Peierls relief (mostly in metals). 
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4 Conclusions 



Perturbation approach allows obtaining linear equation of motion for kink in 
pinning potential. This equation was used for exploring fluctuation-assisted 
depinning of the kink affected by dc and ac forces. The obtained behavior 
of depinning probability on ac force amplitude is rather intriguing. The 
fact that this force can result either in decrease or increase of the probability 
(depending on amplitude of the force and strength of the stopper) means that 
recombination enhancement, modeled in our work by the harmonic ac force, 
can cause both softening and hardening of the crystal. It should be noted that 
both these effects were observed in ionic crystals with low Peierls potential, 
where dislocation segments move (under sufficiently high load) without kink 
formation pp. However, in crystals with covalent and mixed covalent-ionic 
bonds, where the Peierls mechanism of dislocation motion is prevalent, the 
softening is more common phenomena (though the hardening was observed 
too, for example, in GaAs crystals during investigation of photoplastic effect 
[TB]). Our analysis shows, that hardening of Peierls- type crystal due to ac 
excitation can occur in the case when the depinning of migrating kinks is the 
bottleneck of the process. 

It should be noted also that by introducing the ac force we simulated 
only the oscillations of the kink, while the pinning centre was assumed to 
be immobile. Obviously, it is quite possible that the stopper can also be 
influenced by ac excitation, particularly when this center is involved in ther- 
mal recombination process and, consequently, the heat release stimulates its 
oscillations. This can essentially affect the depinning probability |7J. Exten- 
sion of our study to the case when both kink and stopper are subjected to 
ac forces can be the issue of further development. 
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